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OS. 

D ' Abstract. By solving the Fokker-Planck equation directly we examine effects of annihilation, particle escape and 
(J^ I injection on the form of a steady-state positron distribution in thermal hydrogen plasmas with kT < rrieC^. The 
C^ • positron fraction considered is small enough, so it does not affect the electron distribution which remains Maxwellian. 
T-H I We show that the escape of positrons in the form of electron-positron pairs and/or pair plasma, e.g. due to the 
diffusion or radiation pressure, has an effect on the positron distribution causing, in some cases, a strong deviation 
from a Maxwellian. Meanwhile, the distortion of the positron spectrum due to only annihilation is not higher than a 
few percent and the annihilation line shape corresponds to that of thermal plasmas. Additionally, we present accurate 
r — ' formulas in the form of a simple expression or a one-fold integral for energy exchange rates, and losses due to M0ller 
^^ ' and Bhabha scattering, e~^e~-, ee- and ep-bremsstrahlung in thermal plasmas as well as due to Compton scattering 
CN ] in the Klein-Nishina regime. 

^^ ■ Suggesting that annihilation features observed by SIGMA telescope from Nova Muscae and the IE 1740.7-2942 are 
(*^ I due to the positron/electron slowing down and annihilation in thermal plasma, the electron number density and the 
■^ ' size of the emitting regions have been estimated. We show that in the case of Nova Muscae the observed radiation is 
'"^ " coming from a pair plasma stream (ne+ « Ug- ) rather than from a gas cloud. We argue that two models are probably 
I I relevant to the IE 1740.7-2942 source: annihilation in (hydrogen) plasma ne+ ^ n^- at rest, and annihilation in the 
O ' pair plasma stream, which involves matter from the source environment. 

^ \ Key words: diffusion - plasmas - radiation mechanisms: miscellaneous - stars: individual: Nova Muscae 1991 - stars: 
!• ' individual: IE 1740.7-2942 - ISM: general - Galaxy: center - gamma rays: theory 

f 

1. Introduction 

Positron production and annihilation are widespread processes in nature. Gamma-ray spectra of many astrophysical 
sources exhibit an annihilation feature, while their continuum indicates presence of mid-relativistic thermal plasmas 
with kT <i 200 — 300 keV. Spectra of 7-ray bursts and Crab pulsar show emission features in the vicinity of 400-500 
keV (Mazets et al. 1982; Parlier et al. 1990), which are generally believed to be red-shifted annihilation lines. Recent 
observations with SIGMA telescope have exhibited annihilation features in the vicinity of ~ 500 keV in spectra of 
two Galactic black hole candidates, IE 1740.7-2942 (Bouchet et al. 1991; Sunyaev et al. 1991; Churazov et al. 1993; 
Cordier et al. 1993) and Nova Muscae (Goldwurm et al. 1992; Sunyaev et al. 1992). A narrow annihilation line has 
been observed from solar flares (Murphy et al. 1990) and from the direction of the Galactic center (Leventhal et al. 
1978). 

The region of the Galactic center contains several sources which demonstrate their activity at various wavelengths 
and particularly above several hundred keV (e.g., see Churazov et al. 1994). Escape of positrons from such a source or 
several sources into the interstellar medium, where they slow down and annihilate, can account for the 511 keV narrow 
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line observed from this direction. The IE 1740.7-2942 object has been proposed as the most hkely candidate to be 
responsible for this variable source of positrons (Ramaty ct al. 1992). This would only require that a small fraction 
of e+e~-pairs, which is generally believed to be produced in the hot inner region of an accretion disc, escapes into 
surrounding space (Meirelles & Liang 1993). Nova Muscae shows a spectrum which is consistent with Comptonization 
by a thermal plasma kT^ ^ 100 keV in its hard X-ray part, while a relatively narrow annihilation line observed by 
SIGMA during the X-ray flare on 20-21 January, 1991 implies that positrons annihilate in a much colder medium 
(Gilfanov et al. 1991; Goldwurm et al. 1992). 

Numerous studies of positron propagation and annihilation in cold interstellar gas (e.g., see Guessoum et al. 1991 
and references therein) have been inspired by observations of a narrow 511 keV line emission from the Galactic center 
region. Relativistic pair plasmas have been a matter of investigation during a decade (for a review, see Svensson 1990). 
In all thermal models, however, particles are assumed to be Maxwellian a priori and very often one only pays attention 
to the relevant relaxation time scales. Herewith, the annihilation line shape, the main feature which can be actually 
measured, is strongly influenced by the real particle distribution. The latter can differ from a Maxwellian under certain 
circumstances, such as particle injection, escape, and annihilation. It is thus of astrophysical interest to study particle 
distributions in various physical conditions. 

In this paper we use a Fokker-Planck approach to examine the effects of annihilation, particle escape and injection 
on the form of a steady-state positron distribution in thermal hydrogen plasmas with kT < rrieC^. Pairs are assumed 
to be produced in the bulk of the plasma due to 77-, 7-particle, or particle-particle interactions, or to be permanently 
injected into the plasma volume by an external source. We don't touch here upon the cause of particle escape, it 
could be of diffusive origin or due to the radiation pressure (e.g., see Kovner 1984). Since the plasma cloud serves as 
a thermostat, it is therefore reasonable, as the first step into the problem, to consider that the electron distribution 
approaches Maxwellian. The positron fraction considered is small enough, so it does not affect the electron distribution. 

Suggesting that the features observed by SIGMA in ^ 300 keV region are due to the electron-positron annihilation 
in thermal plasma, we apply the obtained results to Nova Muscae and the IE 1740.7-2942 source in order to get the 
parameters of the emitting regions where the annihilation features have been observed. 

In Sect. 2 the Fokker-Planck treatment is considered and we present a method to obtain a steady-state solution. 
The reaction rate formalism is introduced in Sect. 3. The expressions for energy changes and losses due to Coulomb 
scattering, bremsstrahlung and Comptonization are given in Sect. 4-6. Electron-positron annihilation is considered 
in Sect. 7. The results of calculation are discussed in Sect. 8. In the last section (Sect. 9) we discuss the physical 
parameters of the emitting regions in Nova Muscae and the IE 1740.7-2942 source. Throughout the paper units 
h — c — me — 1 are used. 

2. The Fokker-Planck Equation: Positron Spectrum 

Assuming the isotropy of the positron energy distribution function 7(7), the steady-state Fokker-Planck equation takes 
the form 

if {^[^(^)/(t)] - ^(7)/(7)| - [^(7) + i?(7)]/(7) + ^(7) = 0, (1) 

where J f{'j)d'y = 1, 7 is the positron Lorentz factor, P(7) = d-j/dt is the dynamical friction (energy loss rate), 
Z)(7) = d{A'y)'^ /dt is the energy dispersion rate, ^(7) and £^(7) are the annihilation and the particle escape rates, 
respectively, and £'(7) is the positron injection term. 

In the steady-state regime, without sources and sinks, the kinetic coefficients obey the equation (Lifshitz & 
Pitaevskii 1979) which results from the absence of the flux density in the energy space, 

D'il)fi{l) + D{^)I[{l) = ^(7)/i(7), (2) 

where /i(7) is to be a Maxwellian distribution /i(7) ^ 7(7^ — lY''^e'~^'^'^ {kT is the dimensionless plasma tempera- 
ture). This equation fixes a relation between the coefficients 

D{l) = j^J'dih{^')P{^'). (3) 

We emphasize that the coefficients of the Fokker-Planck equation have an additive property. They represent the 
sum of coefficients for various processes which have to be evaluated separately. 
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Although the plasma cloud serves as a thermostat with true Maxwellian distribution, annihilation and sinks distort 
the distribution /i(7). We are thus looking for the solution of eq. (1) in the form fi"/) = fi{"f)g{'j), which gives an 
equation for the unknown function 5(7) (Moskalenko 1995) 

Eliminating .9(7) in favour of ,f{^) yields the integro-differential equation for the distorted function 

/'(7)-/(7){^ + ^-^}^^{/V + i^)/ci7'-f(A + i.)/.7'x/^^.7'}, (5) 

while [Df]' — Pf = at 7 = 1 was assumed (cf. eq. (2)). The last term in cq. (5) follows simply from conservation of 
the total number of positrons 

[A{l)+E(j)]f{j)dj^ S{j)dj, (6) 

1 Ji 

which is always fulfilled if the source function has the form 5'(7) = <S'(7)x J^ {A + E)f d'y' and J^ Sdj — I. A 
regular singular point 7 = 1 in eq. (5) does not lead to any singularity of the solution, which is Maxwellian-like at 
the low-energy part. Equation (6) gives also an idea of physical meaning of term [-E'(7)/(7)], that is the number of 
positrons with Lorentz factor 7 escaping from the plasma volume per 1 sec. The approach can be easily generalized 
to include inelastic processes, stochastic acceleration etc. 

Equation (4) or (5) can be resolved numerically with an algorithm which reduces it to a first-order differential 
equation. Let fi{"/) is the solution obtained after the i-th iteration, then the equation 

/;+r(7)-/.+i(7){^ + :^-^} = ;^{/'(A + i?)/.d7'-/°°(A + £;)/.d7'x^'^d7'}, (7) 

with the initial condition^ /i+i(l) — allows us to get the next approximation /i4.i(7) of the solution. For eq. (4), 
the condition gi^i{l) = gi{l) could be taken. To start the iteration procedure one can use the Maxwell-Boltzmann 
distribution /i(7), although, in some cases, when a solution of eq. (7) deviates strongly from Maxwellian, that causes 
a deviation in normalization during first iterations. Since a solution of eq. (7) fn{l) multiplied by a constant would 
be also a solution, it has to be normalized in the end of iteration process. This algorithm converges quickly and 
gives a good approximation of the solution already after several iterations. The actual signature of the convergence 
could be an equality /^ [A + E)fi-i dj' = /^ {A + E)fi dj' . The combination of functions (a/i + (1 — a).fi-i), where 
a = const ;$ 1, on the place of fi in the right side allows sometimes to get a convergence faster. 

3. Reaction Rate Formalism 

Below we describe a formalism, which further allows us to calculate the annihilation rate, energy losses and energy 
dispersion rate due to Coulomb scattering, bremsstrahlung, and Comptonization. 

The relativistic reaction rate R for two interacting distributions of particles is given by 

R=—^j—Jdn,Jdn2<j{f3,,/3^)il-l3i-f32)Pr, (8) 

where cr(/3i, /32) is the cross section of a reaction, drii and /3i are correspondingly the differential number density and 
velocity of particles of type i in the laboratory system (LS), (3r is the relative velocity of the particles, the factor 
(1 -|- Si2)~^ corrects for double counting if the interacting particles are identical. 

We consider energetic particles which interact with particles of a thermal gas. Let masses of both types of particles 
be equal (m^ = 1). For isotropic distributions, eq. (8) can be reduced to the triple integral over particle momenta. 
Pi ~ Pili, a-iid the relative angle, cosO = pi • P2/P1P2, 

^ = onTl ^ I dpi^hipi) f dp2^f2{p2) I d(cOs6)-lrPra(^r), (9) 

2(l + Ol2JJo 71 Jo 72 7-1 

^ which follows from suggestion f{pj) = fi{'y)g{'y), where /i(7) is a Maxwellian. 
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where rii is the number density of particles of type i in the LS, fi{pi) are the momentum distribution functions, 
!^ dp^pffi{pi) = 1, 

7, = (1 - /32)-i/2 = ^^^2(1 - /3i/32cos0) (10) 

is the relative Lorentz factor of two colliding particles (invariant). 

Putting the relativistic Maxwell-Boltzmann distribution for the electron gas (pay attention to the normalization) 
together with the monoenergetic distribution for the beamed particles, 

e-71/feT 

f2{p2) = ^5{p2-p), (12) 

P2 
into eq. (9) yields 

/I foo 2 

rf(cos0)/ dpi^7r/3.a(7.)e-^^/'=^, (13) 

where Kj is the j-order modified Bessel function. 

Using eq. (10) to eliminate cos0 in favor of p^ = IrPr and changing variables from pi to 71 one can find 

^(^) = o(^^s ^2'KiMl.T^ 2. r '^P^ r ^^1 -<^{Pr)e-'^'"^, (14) 

2[l + di2)kTK2(l/kT)Yp Jo Jy- 7r 

where Ji = 77r(l ± PPr)- After integrating over 71, the reaction rate can be exhibited in the form (Dermer 1985) 

^(^) - (l + .r2)^aAT)7^/? [ ''^ i^^^^ sinh(77.^/3./fcT)e--/-. (15) 

Another form of the reaction rate for interacting isotropic distributions of particles (eqs. [11], [12]) was found useful 
for some purposes (Dermer 1984) 

where 7c = '^^+^^ is the Lorentz factor of the center-of-mass system (CMS), and 7!^ ~ '' ,'' ^^'•"''■' . 

'" V2(7.+l) y \ J^ Ic V2(7. + l) 

If we are interested in energy losses suffered by the energetic particles in an isotropic gas, it is necessary to weight 
the cross section in eq. (14) or (16) by the average LS energy change per collision (A7). The concrete form for (A7) 
depends on the studied process. Hereafter we will consider the reaction rate and energy losses per one positron in the 
unit volume (n2 = 1), while Ue = ni will denote the electron number density. 

4. Coulomb Collisions 

Speaking about the Coulomb scattering one usually implies the lowest order approximation, which is called M0ller 
scattering when referred to identical particles e e*, and Bhabha scattering when referred to distinct particles e+e". 
The effect of bremsstrahlung in ee-collisions is strictly not separable from that of scattering, however, it is convenient 
and generally accepted to treat them separately. Expressions for Coulomb energy losses and dispersion have been 
obtained by Dermer (1985) and Dermer & Liang (1989). Here we describe briefiy their results for the self-consistency 
of consideration. 

The average LS energy change during a collision is (asterisk denotes CMS variables) 



1 f ,■, ,.d^o-^ 



'^^'^^O^/"''"^?^^^' '"' 
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where {d^<y*(j^^Jdp'*^) is the differential cross section, (fp'* = p'*'^dp'*d{cos-ip'*)d(l)'*, ip'* and (f)'* are the polar and 
azimuthal angles, respectively. The LS energy change expressed in these variables is 

A7==7c(7'* -7*)+/3c7c[(p'*cos?A'* -p*)cosw* - p'* sin ?A'*cos0'* sin w*], (18) 

where f3c is the CMS velocity, 

7* - Vhr + l)/2, 



P* = Vhr - l)/2, 



(19) 



are the Lorentz factor and momentum of a particle in the CMS prior to scattering, 7'* and p'* are those after scattering, 
and oj is a kinematic angle 

coscu* ^{f3*-f3,)/p*p,, 
sinuj* = (3i (32 sin ef-yrPrlcPc- 

Energy losses of a particle due to elastic Coulomb scattering are given by eq. (16) with the cross section weighted 
by (A7). Using azimuthal symmetry of the cross section, Dermer (1985) obtains 

"-' "-e^ I J., o ., /o77 T^T / J./* / J/_„„,;.'*^ " "Coul 



d7r /3r7rv2(7r+T) / d^'* I d{cosip'*) 



dt 2{l + Si2)kTK2{l/kTh^pJi " '-' " v '-^' ' ' J ^ J ' ^ ' d-/'*dicos^ 



X J' d7,e"" ""■^'' |^^(y*_^*) + (P_cos^'*-l)(7-7c7*)|, (21) 

since coscj* = (7 — '^cJ*)/jcPcP* ■ In the case of elastic scattering 7* — 7'* and p* = p'* , that gives 
d-f Ue 



dt K2{l/kT)-i^P 



djr Mre'^'^'-'^^Y 



.*2 , l*kT 



-fp*' + ' sinh(77,/3/3,/fcr) - 7/3p*7* cosh(77,/3/3,/fcr) , (22) 

V2(7r + 1)/ I 



where 



^ = 4^ ^(cosr)(l-cosr)^gf^. (23) 

The value of f/'maa; can be assigned from geometrical consideration: tt for distinct particles and 7r/2 for identical particles. 
The minimum scattering angle ip'^^j^ can be related to the excitation of a plasmon of energy ujp. The correction for 
double counting in the case of identical particles appears now as the above condition for i^'^^x ■ 

Integration of eq. (23) with M0ller (e^e^) and Bhabha (e^e^) scattering cross sections (Jauch & Rohrlich 1976) 
gives 

^M. = ^^2^:^:1 ^^2 {^(27*^ - 1)^ (lnA + lnV2 + i) - (27^-7*^ - ^) ln2 + 1(7^ - 1)^} , (24) 

yBH = ^^2^^f^_ ^^2 {^(27*^ - 1)^ In A - ^(227- + I47- - /?- + 6)} . (25) 



The term In A = In \/yz——^pt^ appearing in eqs. (24)-(25) is the Coulomb logarithm. It is a slowly varying function 
of 7*, and often can be approximated by a constant. In the Born regime for the cold plasma limit, the Coulomb 
logarithm is given by Dormer (1985) lnAe±e± = In ( ™^'^ (1 — -)VtH~T), InAgigq: — lnAg±e± + ln\/2. Where the 
plasma frequency LUp can be obtained from the usual expression by replacing the electron rest mass with an average 
inertia per gas particle {j)kT (Gould 1981), ujp = Airr^.c^ne/ {'^)kT- 

Substitution of the Rutherford cross section yields the cold plasma limit 

5^ = -^ In A. (26) 

dt (3 ^ ' 
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The energy dispersion coefficients d{A'j)'^ /dt can be obtained from cq. (3). Anotlier way is to square eq. (18) and to 
follow the above-described method. For M0ller scattering of an electron by a thermal electron distribution the correct 
form of the coefficient has been obtained by Dermer & Liang (1989) 



dt 
where 



^l/kT 



M0 



2kTK2{l/kT)j^p J, 



d-fr 






x^o(li7^-y(7^+ 7*^/3*^) 



77177* (2Xi - X2) + V2 III*' - -f 



27rrf 



12 = 



7*2(^*2 _ 1)2 

27rrg 
7*2(7*2 _ 1)2 



-(27*2 -1)2 + (1-2 In 2) ( 27* 



die Tc^ '=^ 



^(7*2 
12^' 



(27*2 - l)2(lnA + lnA/2) - (^27*4 - 7*2 _ ^\ + 1(^*2 _ ^) 



I2 
2 



(27) 



5. Bremsstrahlung 

Electron-positron bremsstrahlung is a well-known QED process, but the calculation of its fully differential cross 
section for the photon production is very laborious, the resulting cross section formula is extremely lengthy and it was 
obtained quite recently (Haug 1985a,b). In e+e^-colhsions both particles radiate, and that brings some uncertainties 
in calculation of the particle energy loss, increasing particularly as the positron energy closes in the electron gas 
temperature. The exact energy loss rate can be obtained using the cross section differential in the energy of the 
outgoing positrons, but no expression for this quantity is available. As it will be shown, the bremsstrahlung energy 
loss is small in comparison with Coulomb and Compton scattering losses, and that allows us to approximate it by the 
radiated energy rate. We shall, hereafter, speak about the particle energy loss taking into account the above remark. 
An average energy loss through bremsstrahlung is given by 



(A7) 



1 



(7b (Ir) 



dk* 



dn* k- '^''^^* 



dk*dn* 



(28) 



where [d^ a^/ dk* dVt*) is the bremsstrahlung differential cross section in the CMS, and k is the LS energy of the radiated 
photon. It can be expressed as 



(A7) 



IcQc 



CTbilr 



where Qcmilr) — J dk* k*{d(7^/dk*). For e'^e bremsstrahlung Haug (1985c) gives an approximation 

_ 16 2/2 (1.096- 0.523p* + 0.1436p*2 + 1.365p*3_o.532p*4), 7*^8/5; 

Qcm - g ar^ <^ ^^, j^^^, ^^,j _ ^,^g _ Q^^26 + 1.575^*'^ - 0.7967*-2), 7* ;> 8/5, 



(29) 



(30) 



where a is the fine structure constant, and p* , 7* are the CMS variables given by eq. (19). The same for e e 
bremsstrahlung is (Haug 1975) 



Q 



8ar 



2P 



1- 



Ap* 



> 



^]Hi*+P*: 



T I 37* 
Then, starting from eq. (14) and taking into account eq. (29) we get 



(31) 



dprPrilr - l)^'^Qcm{lr)e-''^^"''^ 



d£l_ _ ne 

dt ~ V2(l + (5i2)i^2(l/fcT)72/3 7o 

x{7/37,/3, cosh(7/37,/3./A:T) - (fcT + 7 + 77,) sinh(7/37,/3,/fcT)}. 



(32) 
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In a hydrogen plasma the moving positron suffers energy losses due to e+e^- and ep-brenisstrahlung. For equal e^ 
and p densities, e^e~ bremsstrahlung gives the dominant contribution to the energy loss in the whole energy range. 
At the high energy limit e"'"e~-bremsstrahlung energy loss becomes equal to that of ep and exactly twice the ee energy 
loss; herewith in the Born approximation e+p and e~p cases are identical (Jauch & Rohrlich, 1976). An expression for 
energy loss due to ep-bremsstrahlung was obtained by Stickforth (1961) 

^ - _ 2 2 r87/3[l - (7 - l)/4 + 0.44935(7 - 1)^ - 0.16577(7 - 1)^], 7 ^^S 2; 

dt~ 3 ' ^1/3-^67 ln(27)- 27 -0.2900], 7^2. ^"^"^^ 

6. Compton Scattering 

The presence of photons in a thermal plasma leads to essential energy losses due to Compton scattering. Thomson 
limit remains a good approximation while the photon energy is ^ 1 (the rest mass of the electron) and the electron 
Lorentz factor is not too high. As the photon energy reaches ^0.1 the difference from the classical limit becomes 
large, the principal effect is to reduce the cross section from its classical value. Numerous X-ray experiments show that 
the actual temperature of plasmas in astrophysical sources (far) exceeds 0.05 and the particle Lorentz factor exceeds 
often few units, that is why we consider the Klein-Nishina cross section. 
The particle energy loss rate due to Compton scattering is given by 

d'y 1 l'°° r''* 

' ' dujf^iuj) / dkk<7KN(k){Aj), (34) 



dt 272/3 io 

where 7, (3 are the LS particle Lorentz factor and speed prior to scattering, lo is the initial photon energy in the LS, 
the background photon distribution f^{uj) is normalized on the photon number density n~^ — j duju"^ fj{ui) or on the 
energy density as Uph — J uj'^ J-y{ijj) duj, k^ = ujj{l ± /3), and an average particle energy change due to the scattering is 

(A7) = J- fdk' f d(cos 0') ^p^ A7. (35) 

c^KN J J dk'd[cosO') 

The Klein-Nishina differential cross section in the positron-rest-system (PRS) is expressed in terms of initial k and 
final k' photon energies (Jauch & Rohrlich 1976), 

k' _ 1 

k ~ l + kil-cos9')' 

where 0' is the photon scattering angle in this system. The particle energy change in the LS due to the recoil effect is 

A7 = w-fc'7(l+/3cosp'cos6''), (37) 

where p' is the angle between the incoming photon and positron velocity vectors in the PRS, /3cosp' = (cu/jk) — 1 . 
After the integration one can obtain 

du;f^{u;)[S{j,u;,k+)-S{j,u;,k-)], (38) 



31 5 3 \ , , , , 11, 3 1 1 , ,, 

Y + fc + 2fc^ j '^^'' + ^^ - y ^ - fc + WkTT) + i2(2fc + 1)2 + ^-^(-^^) 

^ + 6 + I) ln(2. + 1) - ^fc + ^^^ - j^^, + 2L^.i-2k)] , (39) 
and Li2 is the dilogarithm 

U (-2k) r\n(l x)^ /E.^l(-2fc)V^^ k < 0.2; 

'^ ^ io ^ a; ~ 1-1.6449341 +iln2(2fc+l) -ln(2fc + l)ln(2fc) + E"i*"'(2fc+l)-% fc> 0.2. 



dj 

dt 


■nrl 

272/3 


L 


where 






Sil,uj 


,fc) = 


-- LO 
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Formulas (38)-(39) give exactly the same result as Jones' (1965) cq. (13). The delta-function approximation of 
the photon distribution /^(a;) ^ S{lj — ujo)/co'^ can sometimes be used for evaluation of the integral (38). We have 
found that in some cases it shows a good agreement with exact calculations, e.g. for the Planck's distribution with 
LOo ^ 2.7kT (see Fig. 2). 

The Thomson limit of the Compton scattering can be obtained similarly by equating fc = fc' in eq. (36) 

S) = -?-e^t/p.7^/3^ (40) 



dtjj, 9 



For the energy dispersion rate one can get 



d(A7)2 



^^Xrl{u^^)j'P'{6^'P^ + l), (41) 

rp 45 



dt 
where (o;^) = j uj'^f^{uj) dio. 

7. Annihilation Rate and Spectrum 

The annihilation rate A (7) for monocnergetic positrons in Maxwell-Boltzmann electron gas can be directly obtained 
from eq. (15) by substitution of the annihilation cross section (Jauch & Rohrlich 1976) 

7^7.2 f I ^ 

""- (>) ^ R2( \u (7r + 4 + -) ln(7, + v/7r' - 1) - Mir + 3) . (42) 

7rPr(7r + 1) I 7r J 

The spectrum of photons dV /dv, which are emitted in the annihilation is given by 

^ = ^+n-|./+(7+)rf7+|/-(7-)^('^,7+,7-)rf7-, (43) 

where v is the dimensionless photon energy, /±(7±) are the arbitrary isotropic particle distributions / f±d'y± = 1, n± 
and 7± are the e* number densities and Lorentz factors. An analytical expression for the angle-averaged emissivity 
per pair of particles, 

r ^*2/3* d<7 

iJ(;.,7+,7-)= /dcosr^-^ — (;.,7+,7_,cosr), (44) 

J "f+-f^ dv 

was obtained by Svensson (1982), here ^ is the differential cross section for emission of a photon with LS energy v. 

8. Calculations and Analysis 

The rates obtained in the paper were integrated over the Maxwellian distribution in order to compare with well-known 
results for the thermal plasma. The annihilation rate was tested with annihilation rate of an e+e^ plasma (Ramaty & 
Meszaros 1981), bremsstrahlung energy losses were compared with e"'"e~-, ee-, and ep-bremsstrahlung luminosities of 
thermal plasmas (Haug 1985c). Two more tests on Coulomb energy losses and bremsstrahlung were carried out with 
calculations by Dermer & Liang (1989). An excellent agreement was found. Compton energy loss eq. (38)-(39) coincides 
with the Thomson limit as a; ^- 0. Besides, we have found that the formulas obtained can be also successfully applied 
for the calculation of the bremsstrahlung luminosity and annihilation rate of the thermal plasma by replacing the 
positron Lorentz factor with the average one over the Maxwellian distribution {^i)kT = 'ikT + [Ki{l/kT)/ K2{l/kT)\. 

The relevant energy loss rates {—d'y/dt) and annihilation rate per one positron are shown in Fig. 1 and 2. All 
values are provided dimensionless, in units UeTrr'^, the Coulomb logarithm was taken a constant In A = 20. M0ller and 
Bhabha energy losses show negligible difference and dominate over the others except Compton scattering, which is quite 
effective and can prevail at large Lorentz factors of positrons (electrons). Low energy particles gain energy in Coulomb 
scattering with thermal electrons that appears as the sign change of d"f/dt. Energy losses due to bremsstrahlung are 
negligible in comparison with others. Annihilation rate is small in comparison with the relaxation rate, so that most 
of positrons annihilate after their distribution approaches the steady-state one. 

The energy losses due to Compton scattering (Fig. 2) have been calculated in the Thomson limit eq. (40) and in the 
Klein-Nishina regime for a Planckian spectrum f-yiuj) = [e^l^^ — 1)^1/2.404(^7)"^, and the i5-function approximation. 
The energy loss rates due to the Comptonization on Planck's photons arc shown for two photon temperatures, the 
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Fig. 1. Shown are the calculated annihilation rate (A), energy losses due to bremsstrahlung (e"'"e~, ee, and ep) as well as 
Coulomb energy losses (C) and dispersion coefficients (D) in thermal hydrogen plasmas. All values are provided dimensionless, 
in units Ueivr^. Low energy particles gain energy in Coulomb scattering with plasma particles that appears as a sign change 
and shown by bold dotted lines. 
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Fig. 2. Energy losses due to the Comptonization. The losses are shown in the units neirr^, the photon number density has 
been taken equal to that of the plasma electrons n-, — Ue- The thin lines (P) show the energy losses on Planck's photons. The 
delta-function approximation of Planck's distribution with ujq = 2.7kT is shown by the dotted lines. Thick solid lines (T) show 
the Thomson limit of the Compton scattering. The dispersion coefficient in the Thomson limit is shown by the solid line (D). 
The inset shows the enlarged low-energy part without the dispersion. 
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Fig. 3. The distorted positron distribution f('y) for an electron temperature kT = 0.1 with and without positron escape 
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with 70=2 (thick lines) and 4 (thin lines). A Maxwell-Boltzmann distribution is shown by a dotted line. The right panel shows 
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Fig. 4. The spectra of photons from annihilation of positrons with Maxwellian electrons for kT = 0.01 and 0.1 with and without 
positron escape i?(7) = 0, 10, 100. The source function of positrons was taken a Gaussian ~ exp[— (7 — 4)^]. The spectra are 
provided dimensionless, in units n-n+yrr^, where n± is the e number density and r^ is the classical electron radius. 



(5-function approximation of Planck's distribution with djg = l.lkT gives similar results. For the clear comparison with 



Fig. 1 the photon number density have been taken equal to that of the plasma electrons n. 



easily generalizing 



for an arbitrary n~^ by trivial vertical shift of the curves. For the coherence, in all calculations the energy density 
of photons C/p/i was taken equal to that of Planck's distribution k, l.lkTn-^. Shown also is the dispersion coefficient 
calculated in the Thomson limit eq. (41). The radiation can provide some heating for the cold particles similar to that 
in the Coulomb scattering. Very low-energy particles gain energy due to Comptonization that appears as a sign change 
of the energy losses (see the inset in Fig. 2). Clearly, the effect results from using the Klein-Nishina cross section. 

At small positron Lorentz factors, the Coulomb energy losses dominate the losses due to Comptonization over 
the variety of photon temperatures and densities (cf. Fig. 1 and 2). Qualitatively it means that high photon density 
leads to the cooling of plasma preferentially through high-energy particles. Herewith, the Coulomb scattering mixes 
particles so that the plasma remains nearly Maxwellian. Therefore the energy losses due to Comptonization would be 
only important for the high-energy tail of the particle distribution, which becomes narrower. The precise shape of the 
distribution would be driven by the balance of income and outcome energy fluxes. 
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Fig. 5. The distorted positron distribution 7(7) for values of electron temperature fcT = 0.1, 0.3, and 0.5 (left panel). The 
positron escape rate was taken the same E = 50/3 for all three cases. Corresponding Maxwell-Boltzmann distributions are 
shown by dotted lines. The right panel shows the spectra of photons from annihilation of positrons with Maxwellian electrons 
for fcT = 0.1, 0.3, and 0.5. The line shapes for annihilation of Maxwellian positrons with Maxwellian electrons are shown by 
dotted lines. 



Positrons could be injected into the hydrogen plasma volume by an external source or produced in the bulk of 
the plasma. In the latter case the form of the source function is governed by the nature of the processes involved. 
Electron-positron pair production in ep-collisions becomes possible when the electron interacting with a stationary 
proton has the Lorentz factor exceeding 3, for ee-collisions one should exceed 7 when one interacting particle is at rest. 
If the pair is to be produced in two-photon collisions, the photon energies, w^, and the relative angle, 0, must satisfy 
the condition lo\lo2 > 2/(1 — cosO). Low plasma temperature is consistent with a small positron fraction in the plasma 
since the positrons could be produced by the relatively small number of head on collisions of energetic photons and/or 
electrons from the tail of Maxwellian distribution. 

If the particle production is not balanced by annihilation it could lead to escape of e^e"-plasma, since the grav- 
itation near a compact object can't prevent pairs from escaping. Two independent mechanisms, at least, diffusion 
and the radiation pressure result in escaping of particles from the plasma volume. We, therefore, explore these factors 
separately. If particles escape due to the radiation pressure, it is natural to suppose that the escape probability E{'-f) 
is a weak function of the particle Lorentz factor, we thus put it a constant. In the case the escape is of diffusive origin, 
the diffusion coefficient is a function of particle speed I? ~ /3. We thus consider two functional forms for the escape 
probability E ^ (3, and E ~ ^/]3 which simulates the case when both mechanisms operate simultaneously. 

Calculations of the distorted function 7(7) have been made (Fig. 3) for the source function in the form of monoen- 
ergetic distribution S{j) — (5(7 — 70), power-law S{'y) — 2/-f^, and Gaussian £'(7) = exp[— (7 — 4:)^]/y/n. The escape 
rate was taken energy-independent E = 0, 10, and 100 in units neTTTg, which is negligible, medium and very high in 
comparison with the time scale of the Coulomb energy losses (cf. Fig. 1). It demonstrates an effect of blowing away of 
(unbound) electron-positron pairs by radiation pressure. 

The behavior of the solution /{"/) depending on the injection function and escape rate is quite clear from the Figure. 
One can show that the right side of the eq. (4) and (5) is negligible at 7 -^ 1, the solution is therefore Maxwellian- like. 
Beginning from some point, the term J. {A + E)f d-y' becomes non- negligible that leads to some increasing of the 
derivative g'{"f) and deviation of the solution from Maxwellian. Thus, a bump is forming. At some Lorentz factor the 
last term in the right side of eq. (4) and (5) is switching on, which leads to some decreasing of the derivative or could 
even change it to a negative value. At large Lorentz factors the right side of the equations again approaches zero (see 
eq. [6]). Generally, if the energy of injected particles essentially exceeds the average one of plasma particles it leads to 
an extended tail, while the correct normalization of the whole solution thus requires some deficit at low energies. 

Typical spectra of photons from annihilation of these positrons with Maxwellian electrons are shown in Fig. 4 for 
electron temperatures kT — 0.01, and 0.1. It is seen that as plasma temperature grows the annihilation line widens, 
its height decreases and distortions of its shape become relatively more intensive. 

Another case is shown in Fig. 5. The distorted functions were calculated for electron temperatures kT = 0.1, 0.3, 
and 0.5 while the escape probability in all cases was taken the same E ~ 50/3 (in units Ueirr'^). The actual values 
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of the escape rate in these cases could be inferred from the value of the integral J^ {A + E)f d-y', which is equal to 
PS 24, PS 35, and « 41, correspondingly. Particle injection was taken monoenergetic with energy equal to the average 
energy of plasma electrons. In all cases, the escape leads to some deficit of energetic particles in the tail of distribution, 
while the particle injection appears as a bump. Although the distributions of positrons differ from Maxwellians, their 
annihilation with thermal electrons does not lead to large distortions of the annihilation line form. This latter is very 
similar to the line from annihilation of two Maxwellian distributions. 

Although only few cases have been discussed, the performed calculations have shown that the functional dependence 
of the escape rate is not very important. In all three cases E = const, ~ /3, and ~ ^/P we obtained similar results for 
the same injection function, the difference appears only at very low temperatures kT ^ 0.05. It is quite clear, since f3 
increases from to ?s 1 in a narrow region 7=1 — 1.4 remaining further a constant. The particle distribution actually 
depends on the value J^{A + E)f d'y, energy of the injected particles and their distribution (cf. Figures 3 and 5). In 
absence of the particle injection, the escape of particles operates as an additional mechanism for the plasma cooling. 

9. Nova Muscae and IE 1740.7 2942 

Recent observations with SIGMA telescope have revealed annihilation features in the vicinity of ^ 500 keV in spectra 
of two Galactic black hole candidates, IE 1740.7-2942 (hereafter the IE source; Bouchet et al. 1991; Sunyaev et al. 
1991; Churazov et al. 1993; Cordier et al. 1993), and Nova Muscae (Sunyaev et al. 1992; Goldwurm et al. 1992). 
During all periods of observation the hard X-ray emission, 35-300 keV, was found to be consistent with the same law. 
Observations of Nova Muscae after the X-ray flare (January 9, 1991) are well fitted by a power law of index 2.4 — 2.5 
or by Sunyaev-Titarchuk (1980) model with kT w 55 — 75 keV and r « 0.4 — 0.5 in the disc geometry, the spectrum 
of the IE source is well described by Sunyaev-Titarchuk model with kT « 35 — 60 keV and t « 1.1 — 1.9. Meanwhile 
soft 7-ray emission of these sources seems to be highly variable. 

During the last 13 hr of a 21 hr observation on January 20-21, 1991, a clear emission feature around 500 keV was 
found in the spectrum of Nova Muscae (Fig. 6), with a line flux of « 6 x 10~^ photons cm~^ s~^, and an intrinsic 
line width ^ 58 keV (Sunyaev et al. 1992; Goldwurm et al. 1992). Since the first 8 hr of the observation did not give 
a positive detection, the inferred rise time is equal to several hours. The next observation, held on February 1-2, did 
not show this feature restricting the lifetime to ^ 10 days. 

The Galactic center region was intensively monitored by SIGMA telescope since 1990. Three times during these 
years a broad excess was observed in the 200-500 keV region (Fig. 6). 

In an observation performed between 1990 October 13 and 14, a spectacular unexpected feature was found in the 
IE emission spectrum, the corresponding flux was estimated at (0.9— 1.3) x 10^^ photons cm~^ s~^, with a line width 
of 180 — 240 keV (Bouchet et al. 1991; Sunyaev et al. 1991). The observations of this region performed two days before 
(on October 10-11), and a few hours after (October 14-15) did not exhibit any spectral feature beyond 200 keV. The 
total duration of this state is estimated between 18 and about 70 hr. 

Seven October 1991 observations have shown an evident excess at high energies, while the source was in a low 
state (Churazov et al. 1993). The excess was observed during 19 days and was not so intensive as in October 1990, 
the average flux was (1.9 ± 0.6) x 10^"^ photons cm~^ s~^ in the 300-600 keV region. 

The 1992 September 19-20 observational session (Sep. 19.42-20.58) showed a feature beyond 200 keV (Cordier et 
al. 1993), which resembles that of 1990 October 13-14. The fine flux was estimated as 4.28^^50 x 10"^ photons cm"^ 
s^^. The previous (Sep. 18.59-19.30) and the next (Sep. 22.57-23.14) sessions did not show any evidence for emission 
in excess of 200 keV, restricting the lifetime of the state between 27 and about 75 hr, while the rise time approaches 
probably several hours. 

The spectral features observed by SIGMA are, commonly believed, related to electron-positron annihilation. Rel- 
atively small line widths imply that the temperature of the emitting region is quite low, kT w 35 — 45 keV for IE 
and 4-5 keV for Nova Muscae. Since the hard X-ray spectra < 300 keV showed no changes, most probably that 
electron-positron pairs produced somewhere close to the central object were injected into surrounding space where 
they cool and annihilate. Radiation pressure of a near-Eddington source alone can accelerate e+e^-plasma up to the 
bulk Lorentz factor of 70 ^ 2 — 5 (Kovner 1984), while Comptonization by the emergent radiation field (Levich & 
Syunyaev 1971) could provide a mechanism for cooling the pairs which further annihilate "in flight" (for a discussion 
see also Gilfanov et al. [1991, 1994]). If there is enough matter around a source, then particles slow down due to 
Coulomb energy losses and annihilate in the medium. We explore further this last possibility by checking whether the 
inferred parameters of the emitting region are consistent with those obtained by other ways. We assume single and 
short particle ejection on a timescale of hours. It seems reasonable: since the ejection would probably impact on the 
whole spectrum, longer spectral changes would be observable. 
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Fig. 6. Energy spectra of tlie IE source (Boucliet et al. 1991; Churazov et al. 1993, 1994; Cordier et al. 1993) and Nova Muscae 
(Goldwurm et al. 1992) observed by SIGMA are shown together with fits of the authors. For September 1992 flare shown 
is counts s~^ lieV~^. The dashed line in the upper left panel shows the annihilation line shape for Gaussian- like injection, 
~ exp[— (7 — 4)^], of energetic particles into the thermal plasma of kT = 35 keV for E/A = 20. The line is shifted left to 
approach the data. 



Suggesting that the energetic particles slow down due to Coulomb scattering in the surrounding matter, one can 
estimate its (electron) number density 



70 



;A,: 



(45) 



where 70 is the initial Lorentz factor of the plasma stream, c is the light speed, and A^ is the characteristic time scale 
of the annihilation line appearance. The Coulomb energy loss rate in a medium of kT < 0.1 is (dj/dt) w 70 — 100 (see 
Fig. 1). Taking a reasonable value for the bulk Lorentz factor, 70 ~ 3 (e.g., Kovner 1984), one can obtain estimations of 



the order of magnitude as n_ 
for Nova Muscae. 



(Ai/2 days) ^ for the IE source, and n_ w 1.5 x 10* cm ■^ (Ai/5 hr) 
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Table 1. Observational data and parameters of the emitting region. 



IE 1740.7-2942 



1990 Oct. 13-14 1992 Sep. 19-20 



Nova Muscae 



Annihilation rise time, Ai ^2 days (1-2 hr)* few hours ~ 5 hr 

Annihilation lifetime, Ad 18-70 hr 27-75 hr <y 10 days 

Annihilation photon flux, Fsoo (photons cm"^ s"^) 10"^ 4.3 x 10"^ 6 x 10"^ 

Total line flux, L500 (photons s"^) 8.6 x 10"*^ (at 8.5 kpc) 3.7 x 10*^ 7.2 x lO"^ (at 1 kpc) 

Line width, W (keV) 240 180 40 

Plasma temperature, kT^ (keV) 35 — 45 3 — 4 

Column density, Nh (cm'^) ~ 10^^ ~ 10^^ 

Coulomb energy loss rate, d-y/dt 70 100 

Annihilation rate, A 11 

Electron number density, n_ (cm^^) (5 - 20) x 10** 1.5 x 10* 

Size of the emitting region, A (cm) (1.1 - 20) x 10^* (1.3 - 50) x lO" 

* Our estimation. 

If the energetic particles were injected into the medium only once, then the annihilation feature lifetime A^ is 
directly connected with annihilation rate as A^ w Trr'^cn^A{'y). It yields one more estimation of the number density 
in the emitting region 



1.55 X 10^ cm"-'' -— ^ ) . (46) 



TTr^c Aj^ A{"/) \ 1 day 

Annihilation rate ^(7) is a weak function of 7 (see Fig. 1) and we can take it equal to a constant A — A{1) « 1. 
Total duration of the hard state is A^ w 18 — 70 hr for the IE source and Ad < 10 days for Nova Muscae, that 
gives n_ « (5 — 20) x 10* cm"^ and n_ w 1.5 x 10* cm"^ (A^/IO days)"^, correspondingly. The values obtained from 
eq. (45)-(46) restrict the electron number density in the volume where particles slow down and annihilate. 
Being equated eq. (45)-(46) give an obvious relation between the time scales 

^ = ' ^. (47) 

A, A(7o-l)dt ^ ' 

Therefore, to be consistent with the annihilation lifetime the annihilation rise time for the IE source should be 
Aj « 1 — 2 hr. This is supported by the 1992 September 19-20 observation when the annihilation rise time was 
restricted by a few hours. 

The size of the emitting region A could be estimated from a simple relation n^X^ ^ AdL^Qof^ if we assume the 
upper limit for the positron number density n^ < n_. It gives A <; 1.34 x 10^^ cm (Ajj/l day)^/'' w (1.1 — 2.7) x 10^^ 
cm for IE and A <; 1.3 x 10^^ cm (A^/IO days)^/'' for Nova Muscae^, which are well inside of the upper limits 
A < cAi w 2.2 X 10^^ cm (Ai/2 hr) and < 5 x 10^^ cm, correspondingly. From the above consideration follows that 
emitting regions in both sources are optically thin and do not affect the Comptonized spectra at < 300 keV nor the 
annihilation line form. Experimental data and the estimated parameters are summarized in Table 1. 

The column density of the medium where injected particles slow down and annihilate should exceed the value 
Nh ^ An_, which follows from previous estimations for A and n_, viz. 2.1 x 10^^ cm~^ (^d/l day)""'^/^ ^ Nh < 
cAiTi- w 1.1 X 10^^ cm"^ for IE, where we took into account eq. (47), and Nh <; 2 x 10^^ cm~^ (Ad/10 days)"^/^ 
for Nova Muscae. The total column density of the gas cloud measured along the line of sight, where the IE source 
embedded, is high enough Nh w 3 x 10^^ cm"^ (Bally & Leventhal 1991; Mirabel et al. 1991). Note that recent ASCA 
measurements of the column density to this source give a best fit value Nh w 8 x 10^^ cm"^ (Sheth et al. 1996). For 
Nova Muscae the corresponding value is Nh ^ 10^^ cm"^ (Greiner et al. 1991), less or marginally close to the obtained 
lower limit. If, on contrary, one suggests n+ <C n_ it yields a condition Nh 3> 2 x 10^^ cm~^ (A^/IO days)"^/"', which 
considerably exceeds the measured value. 

These estimations put us on to an idea that the 500 keV emission observed from Nova A4uscae was coming from 
e^e^-plasma jet (n-|_ ss n_) rather than from particles injected in a gas cloud'' (n-|_ <C n_), therefore, particles have 



^ We took n+ < n- « 1.5 x 10* cm"^ (Ad/lO days)-^ 

A possibility that Nova Muscae lies in front 
injected into this cloud, away from the observer. 



^ A possibility that Nova Muscae lies in front of a large gas cloud can not be totally excluded. In this case, particles could be 
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to annihilate "in flight" producing a relatively narrow line blue- or red-shifted dependently on the jet orientation. 
If so, then our estimation of the electron number density n_ from annihilation time scale is related to the average 
electron/positron number density in the jet, its total volume is of A*^ ~ 2 x 10'^^ cm"^ (A^/IO days)^. The reported 
6%-7% redshift of the line centroid (Goldwurm et al. 1992; Sunyaev et al. 1992) supports probably the annihilation- 
in-jet hypothesis, although authors noted that statistical significance of this shift is not very high. The large size of the 
emitting region and a small width of the line, both except the gravitational origin of the redshift, since in this case the 
annihilation region have to be quite close to the central object ~ lORg where typical flow velocities should result in a 
much broader line (Gilfanov et al. 1991). The Compton scattering of the anisotropic emergent radiation could provide 
effective mechanism for blowing away and acceleration of e+e^-pair plasma (e.g., Kovner 1984; Misra & Melia 1993) 
cooling it at the same time. Since the maximal energy during the X-ray flare of Nova Muscae released near ~ 1 keV 
(Greiner et al. 1991), the average kinetic energy per particle should be nearly the same (which is consistent with the 
small line width). 

The case of the IE source is not definitively clear, because our estimations give n_ <; n^ in the emitting region. 
Two flares, October 1990 and September 1992, have shown very similar time scales, spectra and photon fluxes, which 
are consistent with single injection of energetic particles into the thermal (hydrogen) plasma. Meanwhile, the redshift of 
the line ~ 25% reported by authors (Bouchet et al. 1991; Sunyaev et al. 1991; Cordier et al. 1993) implies that positrons 
probably annihilate in a plasma stream moving away from the observer. The estimation of the size of the emitting region 
ruled out its gravitational nature, since it is too large in comparison with gravitational radius of a stellar mass black 
hole. A natural explanation of this controversial picture is that the propagating plasma stream captures matter from 
the source environment and annihilation occurs in a moving plasma volume. In this case the estimation of the electron 
number density n_ is related to the average electron number density in the jet, A^ ^ 2.4 x 10"^^ cm'^ (^d/l day)^ gives 
its total volume, and the jet length has to be of the order of ^ 0.2cAd ~ 5.2 x 10^^ cm (A^/l day). 

While a part of the e"'"e~-pair probably annihilate in a thermal plasma near the IE source producing the broad 
line, the remainder could escape into a molecular cloud, which was found to be associated with the IE source (Bally 
& Leventhal 1991; Mirabel et al. 1991). The time scale for slowing down"* due to the scattering could be obtained 
from eq. (45). Taking ~ 10^ cm~^ for the average number density of the molecular cloud near IE (Bally & Leventhal 
1991; Mirabel et al. 1991) one gets A, ^ 1 year, the same as that obtained by Ramaty et al. (1992). The size of the 
turbulent region in the cloud caused by propagation of a dense jet should be of the same order. It agrees well with 
the length 2-4 ly (15-30 arcsec at the 8.5 kpc distance) of a double-sided radio jet from the IE source found recently 
with the VLA (Mirabel et al. 1992). 

If the lines from the IE source (Fig. 6) were produced by continuous injection of energetic particles, then the 
observations of the narrow 511 keV line emission from the Galactic center allows to put an upper limit on the particle 
escape rate into the interstellar medium. Recent reanalysis of HEAO 3 data has shown that under suggestion of a single 
point source at the Galactic center narrow line intensities are i^sn = (1.25 ±0.18) x 10"'^ photons cm~^ s^^ for the fall 
of 1979 and F511 = (0.99 ±0.18) x 10"^ photons cm'^ s'^ for the spring of 1980 (Mahoney et al. 1994). Taking tq = 1 
yr for the positron lifetime in 10^ cm~^ dense cold molecular cloud (Ramaty et al. 1992), and suggesting one hard state 
of Ad ^ 2 days long per period tq, one can obtain an escape rate E/A « p'"^^'^' ^ 20, where we took F500 = 10^^ 
photons cm^^ s^^ (see Table 1). This is consistent with the upper limits of 1990 October 13-14 spectrum and the two 
most energetic points in 1992 September 19-20 spectrum. The dashed line in 1990 October 13-14 spectrum (Fig. 6) 
shows the annihilation line shape for Gaussian-like injection, ~ exp[— (7 — 4)^], of energetic particles into the thermal 
plasma of kT = 35 keV for E/A = 20. The longest hard state (^ 19 days) with the average flux of F500 ~ 2 x 10^"^ 
photons cm^^ s^^ observed in October 1991 places the upper limit at almost the same level E/A w 10. 

10. Conclusion 

We have presented the accurate formulas in the form of a simple expression or an one-fold integral for the energy 
losses and gains of particles scattered by a Maxwell-Boltzmann plasma. The processes concerned are the Coulomb 
scattering, e^e~-, ee- and ep-bremsstrahlung as well as Comptonization in the Klein-Nishina regime. 

The problem of positron propagation is treated in a Fokker-Planck approach, which can be easily generalized to 
include inelastic processes, stochastic acceleration etc. We have shown that the escape of positrons in the form of 
pair plasma has an effect on the positron distribution causing, in some cases, a strong deviation from a Maxwellian. 
When the energy of injected particles essentially exceeds the average one of plasma particles, the deviation appears 
as a deficit at lower energies and an extended tail of the distribution that leads to a widening and smoothing of the 



* The corresponding annihilation lifetime Ad (eq. [46]-[47]) was obtained for thermal plasma and is not valid for the cold 
medium where positrons mostly annihilate in the bound (positronium) state. 
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annihilation feature in the spectrum. In the ease where the energy of particles injected is close to the average energy of 
plasma particles, the deviation appears as an injection bump and a deficit in the tail of the distribution. Meanwhile, 
it does not lead to visible distortions of the annihilation line shape which is similar to that of thermal plasmas. 

The performed calculations allow us to obtain reliable estimations of the electron number density and the size of 
the emitting regions in Nova Muscae and the IE 1740.7-2942 source, suggesting that spectral features in 300-600 keV 
region observed by SIGMA telescope are due to the electron-positron annihilation in thermal plasma. We conclude 
that in the case of Nova Muscae the observed radiation is coming from a pair plasma jet, n_(- w n_, rather than 
from a gas cloud. The case of IE 1740.7-2942 is not definitively clear, n+ ^ n_. Although the observational data 
are consistent with annihilation in (hydrogen) plasma at rest, the redshift of the line suggests that it could be also a 
stream of pair plasma with matter captured from the source environment. 
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